Linear chain structure of four-a clusters in 16 
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We investigate the linear-chain configurations of four-a clusters in 16 using a Skyrme cranked Hartree- 
Fock method and discuss the relationship between the stability of such states and angular momentum. We 
show the existence of a region of angular momentum (13-18 Ti) where the linear chain configuration is 
stabilized. For the first time we demonstrate that stable exotic states with a large moment of inertia 
(h 2 12® ~ 0.06-0.08 MeV) can exist. 



PACS numbers: 21.60.Jz, 21.30.Fe, 21.60.Cs 

Strong nuclear deformations provide an excellent frame- 
work in which to investigate the fundamental properties of 
quantum many-body systems. Strongly deformed nuclei 
have been identified by the observation of y-ray cascades 
typical of rotational bands. Since the first observation of 
such bands [1], strongly deformed states with an aspect ra- 
tio 1:2 have been found in various nuclei. These bands 
are called superdeformed bands. Furthermore, the hyper- 
deformed bands, in which the deformation is around 1:3, 
have been reported in several experiments [2]. At first, 
those strongly deformed states were found in the heavy nu- 
clei. Such new data have triggered interest in whether more 
exotic states exist in light nuclei where a strong deforma- 
tion above 1:3 could be possible due to a-cluster structure. 

Experimental candidates for strongly deformed states 
with an aspect ratio above 1:3 have been suggested in light 
AN nuclei. One candidate is the four-a linear chain band 
starting around the four-a threshold energy region in lfi O 
suggested by Chevallier et al. [3] in the 12 C (a, 8 J3e) 8 Be 
reaction and was supported by Suzuki et al. Freer et 
al. JH] performed the 12 C( 16 0,4a) reaction and obtained 
a smaller moment of inertia, about 2/3 of Ref. fl. Re- 
cently it has been suggested that these have a loosely cou- 
pled four-a structure Jg], in connection with the gas-like 
0+ (Hoyle) state in 12 C 0]. Another candidate is a six- 
a linear chain state, which has been extensively studied 
both theoretically and experimentally 1181- 41 lh . Wuosmaa 
et al. [flOfl and Rae et al. [11] suggested that the molecular 
resonance state observed in the inelastic reaction 12 C( 12 C, 
12 C(02)) 12 C(02) might be a candidate for the six-a linear 
chain state. Hirabayashi et al. claimed lfl2ll that this has 
a loosely coupled 3a+3a configuration rather than a lin- 
ear chain. The seven-a linear chain state in 28 Si was not 
observed lfl3ll . Despite many efforts no clear experimen- 
tal evidence of a stable a linear chain structure has been 
confirmed and its existence remains an open problem. 

The stability of such linear chain states has often been 
studied through the analysis of small vibrations around 
the equilibrium configuration and with the axial symme- 



try However, it was shown that bending motion 

is an essential path for the transition to low-lying states 
lfl5ll . Thus it is necessary to calculate the stability in a wide 
model space, which contains lower excited states. Two 
mechanisms are important for stabilizing the linear chain 
state. The first mechanism is the quantum-mechanical or- 
thogonality condition to other low-lying states. The sec- 
ond is, as discussed by Wilkinson [ 16], the competition be- 
tween the nuclear attractive and centrifugal forces due to 
rotation of the system: a large moment of inertia such as in 
the linear-chain configuration is favored with a large angu- 
lar momentum. On the other hand, high angular momenta 
would lead to fission of the linear chain due to the strong 
centrifugal force. Detailed investigations are necessary for 
the existence of a region of angular momentum where the 
linear chain configuration is stabilized. 

In this Letter we show that a region of angular momen- 
tum (13-18 h) where the four-a linear chain configuration 
is stabilized exists in 16 0. The cranked Hartree-Fock (HF) 
method is used to investigate the stability of the configu- 
ration and the moment of inertia. The mechanism of the 
stabilization against the decay with respect to bending mo- 
tion and fission and its angular-momentum dependence is 
clarified. 

Until now, most of the theoretical analyses of the linear 
chain structure have been performed using the conventional 
cluster model with effective interactions, whose parameters 
are determined to reproduce the binding energies and scat- 
tering phase shifts of the clusters. Thus it is highly desir- 
able to study the presence of exotic cluster configurations 
based on different approaches, such as mean field models. 
The effective interactions used in the mean field models 
are determined in a completely different way; they are de- 
signed to reproduce various properties of nuclei in a wide 
mass range. The appearance of cluster structure as a re- 
sult of calculations with such interactions and model spaces 
would give more confidence in their presence. Recent de- 
velopments of three-dimensional calculations with Skyrme 
forces enable us to describe both shell-like and cluster-like 



2 




FIG. 1. (Color online) Surface of the total nucleon density distri- 
bution for the initial twisted wave function. We chose the surface 
as the half of the total density. 



configurations, and the interplay between these structures 
based on this approach has been successfully investigated 

fH. 

To investigate the four-a linear chain state in the rota- 
tional frame, we perform cranked HF calculations. We 
self-consistently calculate the cranked HF equation, given 
by 8 (H - u)J) = 0, where H is the total Hamiltonian, to 
is the rotational frequency, and J is the angular momen- 
tum around the y axis. We represent the single-particle 
wave functions on a Cartesian grid with a grid spacing of 
0.8 fm. The grid size is typically 24 3 for ground states 
and 32 x 24 2 for superdeformed states. This accuracy was 
seen to be sufficient to provide converged configurations. 
The numerical procedure is the damped-gradient iteration 
method 11181 . and all derivatives are calculated using the 
Fourier transform method. 

We take three different Skyrme forces which all perform 
very well concerning nuclear bulk properties but differ in 
details: Sly6 as a recent fit which includes information on 
isotopic trends and neutron matter 11911 . and SkI3 as well 
as SkI4 as recent fits which map the relativistic isovec- 
tor structure of the spin-orbit force If2fjll . The force SkI3 
contains a fixed isovector part analogous to the relativis- 
tic mean-field model, whereas SkI4 is adjusted allowing 
free variation of the isovector spin-orbit term. Thus all 
forces differ somewhat in their actual shell structure. Be- 
sides the effective mass, the bulk parameters (equilibrium 
energy and density, incompressibility, and symmetry en- 
ergy) are comparable. 

Here we discuss the stability of the four-a linear chain 
configuration in the rotating frame for ie O. To this end, 
we perform the cranked HF calculations with various ro- 
tational frequencies, oj. For the initial wave function, we 
chose the z-axis as the principal axis and use the twisted 
four-a configuration, as shown in Fig. [TJ We also show 
the corresponding two-dimensional plot in Fig. [2(a). Note 
that this initial is a three-dimensional four a configuration, 
which facilitates the transition of the initial state to low- 



lying states including the ground state during the conver- 
gence process. This was demonstrated for the carbon chain 
states in Refs. 112 ll 12211 . We calculate the rigid-body mo- 
ment of inertia, 0, using the total nucleon density at each 
iteration step. We only consider rotation around the y axis 
(perpendicular to the deformation axis z). 

We first investigate the convergence behavior of the HF 
iterations. To check this, we calculate the coefficient of 
the rotational energy, given by h 2 /2&, at each iteration 
step. Figure [3] shows the calculated results with various 
rotational frequencies versus the iterations in the case of 
the SkI3 interaction. The initial state with the twisted lin- 
ear chain configuration is not the true ground state of the 
HF model space and the solution changes into the true 
ground state after some large number of iterations; how- 
ever the situation depends on the value of the rotational 
frequency oj. In Fig. we see that the rotational frequen- 
cies oj = 0.5, 1.0, and 1.5 MeV/7i (the dashed, dotted, and 
dot-dashed lines, respectively) lead to the true ground state. 
Note that the rigid-body moment of inertia, 0, is a classical 
value, which does not become to be zero even with spher- 
ical shapes. The corresponding density distribution at the 
15000th iteration is plotted in Fig. |2] (b). The frequency 
oj = 0.0 MeV/h (the solid line in Fig. [3) leads to the quasi- 
stable state (see Fig. p2] (c)). At around oj = 2.0 MeV/fi, 
we obtain the state (the thick solid line in Fig. [3J with the 
four-a linear chain configuration, as shown in Fig. [2] (d), 
whereas fission occurs above those rotational frequencies 
(the dot-dot-dashed line in Fig. [3]). 

We next estimate the range of the rotational frequencies 
where the four-a linear-chain configuration can be stabi- 
lized. Figure |4] shows the coefficient of the rotational en- 
ergy, h 2 /2&, versus the rotational frequency oj with various 
Skyrme interactions. We find stable states for the four-a 
linear chain configuration for all of the interactions. For 
the SkI3 interaction, we obtain the lower and upper bounds 
of the rotational frequencies as 1.8 and 2.2 MeV/h. Be- 
tween these the four-a linear chain configuration is stabi- 
lized. The values are 1.9 and 2.2 MeV/h for the SkI4 in- 
teraction and 2.0 and 2. 1 MeY/h for the SLy6 interaction, 
respectively. In these frequency regions where the linear 
chain configuration is stabilized, we can define the rigid- 
body moments of inertia, which are calculated as 0.065 
MeV for the SkI3 and SkI4 interactions and 0.06 MeV 
for the SLy6 interaction. These values are consistent with 
0.063 in a naive picture of rigid-body four-a's laid in linear 
chain with 12 fm as in Fig. 2(d). 

We also estimate the corresponding angular momentum 
where the four-a linear chain configuration is stabilized. 
We calculate the angular momentum using the rigid-body 
moment of inertia obtained and compare it with the value 
calculated by the cranking method. The angular momen- 
tum with the rigid-body moment of inertia, / r id, is calcu- 
lated as 7 r id = &oj. The angular moment calculated us- 
ing the cranking method, 7 cra , is given by 7 cra =< J >, 
where < J > is the expectation value of the angular mo- 
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FIG. 2. (Color online) Total nucleon density distribution calculated using the cranking method for (a) the initial wave function, (b) the 
ground state, (c) the quasi-stable state, and (d) the four-a linear chain state. The isolines correspond to multiples of 0.02 fm~ 3 . We 
normalize the color to the density distribution at the maximum of each plot. 



mentum in the cranking equation. Figure [5] shows the an- 
gular momentum obtained versus the rotational frequency. 
We see that the calculated angular momentum using the 
rigid-body moment of inertia agrees well with that of the 
cranking method, indicating that the rigid-body approxi- 
mation is reasonable for the four-a linear chain states. We 
find that the lower and upper bounds of the angular mo- 
mentum where the four-a linear chain configuration is sta- 
bilized are about 13 and 18 Ti for the SkI3 interaction, 14 
and 18 fi for the SkI4 interaction, and 16 and 18 Ti for the 
SLy6 interaction, respectively. With a such high angular 
momentum, very exotic configuration of the four-a linear 
chain can be stabilized. Fission occurs beyond this angular 
momentum region. Furthermore, it is possible that states 
with even lower angular momenta are stabilized, when the 
coupling effect with low-lying states is taken into account. 
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FIG. 4. Coefficient of the rotational energy h 2 /2& as a function 
of rotational frequency a>. The lines correspond to the different 
Skyrme forces as indicated. 
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FIG. 3. Coefficient of the rotational energy, h 2 /2&, calculated us- 
ing the cranking method versus the HF iterations with various ro- 
tational frequencies w. The symbols (b), (c), and (d) correspond 
to the density distributions given in Fig. [2] 



FIG. 5. Angular momentum as a function of rotational frequency 
w for the Skyrme forces. The lines with solid symbols denote the 
calculated results for the rigid-body moment of inertia, while the 
open symbols denote the results for the cranking method. 
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FIG. 6. Calculated energies of the four a linear-chain states ver- 
sus the angular momentum. The lines correspond to the different 
Skyrme forces as indicated. 



Finally in Fig. [6] we show the energies of calculated four 
a linear chain states versus the angular momentum J. The 
energy is very high, at around -70 MeV with J ~ 14. If 
we extrapolate the energy to lower angular momentum, the 
band head energy at J = is estimated to be around -82 
MeV, which is much higher than the value suggested in 
Ref. H. 

This rather high excitation energy in the Skyrme 
forces is also reported in Ref. [9|]. 

In summary, we have investigated the stability of the 
four-a linear chain configuration in 16 using the Skyrme 
mean field method with cranking. Even if the bending path 
is opened in the three-dimensional space, we obtained re- 
gions of rotational frequency where the linear chain config- 
uration is stabilized. Below this region, the state converges 
to low-lying configurations, with fission occurring beyond 
this region. The frequency range corresponds to angular 
momenta of 13-18 fi. Furthermore, when the coupling ef- 
fect with low-lying states is taken into account there is a 
possibility that states with even lower angular momentum 
are stabilized. The coefficient of the rotation (h 2 /2&) of the 
four-a linear chain configuration obtained is around 0.06- 
0.08 MeV. We have, for first time, shown that states with 
such large moments of inertia are possible in light nuclei 
under conditions of large angular momenta. 

As shown in this Letter, the exotic four-a linear-chain 
state can indeed exist in 16 0. We also investigated whether 
such state can be accessible in the 8 Be + 8 Be fusion reac- 
tion using the time-dependent HF method with the same 
Skyrme force as the present study [23]. We obtained 
a quasi-stable state with a similar moment of inertia, as 
shown here. The HF method is a powerful tool for inves- 
tigating both the static and dynamical properties of nuclei 
in the consistent framework. This method is a promising 
tool to reveal the existence of states with more exotic geo- 
metric configurations, such as longer linear chain 12411 and 



polygons 11611 . which is still an open question for their ex- 
istence. 
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